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We study the influence of boundary conditions transverse to the transport direction for disordered 
mesoscopic conductors both at the Anderson metal-insulator transition and in the metallic regime. 
We show that the boundary conditions strongly influence the conductance distribution exactly at 
the metal-insulator transition and we discuss implications for the standard picture of one-parameter 
scaling. We show in particular that the scaling function that describes the change of conductance 
with system size depends on the boundary conditions from the metallic regime up to the metal- 
insulator transition. An experiment is proposed that might test the correctness of the one-parameter 
scaling theory. 
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I. INTRODUCTION 

More than fourty years after its discovery by Ander- 
son [[jj the disorder-induced metal-insulator transition 
is still the subject of strong theoretical as well as exper- 
imental research ||. One of the major achievements in 
the long history of the Anderson metal-insulator tran- 
sition (MIT) is the renormalization group theory, which 
has also become known as one-parameter scaling theory 
Its basic assumption is that close to the transition 
the change of the dimensionless conductance g with the 
sample size L depends only on the conductance itself and 
not separately on energy, disorder, the size of the sam- 
ple, its shape, the elastic mean free path l D , etc.. Many 
predictions, like the lower critical dimension, or the crit- 
ical behavior 0|| were successfully based on this theory, 
as well as an enormous amount of numerical work that 
aimed at the direct calculation of the scaling function 
(3(g) = dhig/dlnL. Another important consequence of 
the one parameter scaling theory is the prediction of a 
universal conductance distribution P*{g) exactly at the 
MIT Q . Earlier numerical work on the three dimensional 
Anderson model seemed to confirm the universality of 
the conductance distribution JtJ . The dependence on the 
universality class was stressed in ||. 

Recently, however, some doubts have been cast on 
whether the conductance distribution is universal within 
the same universality class. Two different numerical 
studies reported two different forms of P* (g) for the same 
system and it was found that the difference origi- 

nates in the use of different boundary conditions (BCs) 

The idea that P*(g) might depend on the BCs ap- 
pears indeed very natural after the discovery that spec- 
tral statistics, and in particular the energy level spacing 
distribution P(s) exactly at the MIT, do depend on the 



BCs [fL2| . Samples with periodic boundary conditions 
show a much stronger level repulsion than samples with 
hard walls (Dirichlet boundary conditions). 

In this work we show with a numerical analysis of the 
conductance distribution at the critical point that P* (g) 
does indeed depend on the BCs applied perpendicular 
to the transport direction. Choosing the appropriate 
boundary conditions, we can reproduce both the results 
of refs. H and ||. In particular, the average critical 
conductance g c depends on the BCs. This alone already 
implies a dependence of (3(g) on the BCs since g c is de- 
fined as (3(g c ) = 0. We confirm the BC-dependence of 
(3(g) analytically by reinvestigating its form in the metal- 
lic regime with the help of a 1/g expansion. Much to our 
surprise we find that earlier analyses overlooked the effect 
of the BCs by approximating a sum over diffusion modes 
by an integral. Evaluating the sum more carefully, we 
not only find a dependence on the BCs, but also a so 
far unkonwn \n(l e / L) / g term in (3(g) in three dimensions 
that makes (3(g) non-universal in the metallic regime. 



II. NUMERICAL INVESTIGATION AT THE 
ANDERSON TRANSITION 

The model studied is the three dimensional tight bind- 
ing Anderson Hamiltonian with diagonal disorder on a 
simple cubic lattice, 
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The ei are distributed uniformly and independently be- 
tween —w/2 and w/2. The notation < ij > means next 
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nearest neighbors, u is the hopping matrix element which 
we set equal to unity in the following, and w is the disor- 
der parameter. The last sum in eq.(|l]) links correspond- 
ing sites on opposite sides of the cubic sample perpen- 
dicular to the y and z directions, assuming that trans- 
port occurs in the x-direction. Hopping between these 
boundary sites arises when the system is closed to a ring 
(c = 1) and includes a phase factor e l27r< ^, where cj> is the 
magnetic flux in units of h/e inclosed by the ring. Hard 
wall (Dirichlet) BCs correspond to c = 0. The model (|l|) 
shows a MIT at the critical disorder w c ~ 16.5 

The numerical calculation of the conductances uses a 
standard Green's function recursion technique that 
yields the transmission matrix t of the sample. The latter 
is connected to the two-probe conductance of the sample 
by the Landauer-Buttiker formula 



g = tr tt + 



(2) 



where g = G/(e 2 /h) denotes the conductance G in units 
of the inverse of the von Klitzing constant h/e 2 . Whether 
the two-probe conductance formula or the four-probe 
conductance formula is used, is quite irrelevant at the 
metal insulator transition, since the bulk resistance al- 
ways dominates largely over the contact resistance [ p"5| |. 
All conductances were calculated at energy E ~ 0. The 
number of conductances used for each BC and system 
size ranged between 10 5 for L = 6 and L — 8 to 2 • 10 3 
for L = 16. All system sizes L are measured in units of 
the lattice constant. 




FIG. 1. Critical conductance distribution for periodic and 
hard wall boundary conditions, and different sample sizes. 
Sample sizes L = 8, 10, 12, 16 are denoted by circles, squares, 
diamonds and triangles, respectively; open symbols indicate 
periodic BCs, full ones hard walls. The full lines are averages 
over the above system sizes. 

Our main numerical result is shown in Fig.[j], where 
we have plotted the distributions of the logarithm of the 
conductance at the transition for periodic and hard wall 
(HW) BCs, and different system sizes. For the same BC 



the distribution is almost independent of the system size, 
as is to be expected from the criticality of the ensemble 
at w c = 16.5. But the distributions are clearly very dif- 
ferent for the two BCs. The maximum of the distribution 
is considerably more pronounced for periodic BCs than 
for hard walls. A more detailed statistical analysis is pre- 
sented in table | and for the average values (In g) in figure 
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0.377 
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0.392 
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0.334 
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P 


0.413 


0.336 


-1.329 


1.092 
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0.313 


0.306 


-1.777 
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HW 


0.326 


0.310 


-1.710 


1.252 


10 


HW 


0.331 


0.312 


-1.685 


1.246 


12 


HW 


0.338 


0.311 


-1.675 


1.211 


16 


HW 


0.348 


0.319 


-1.614 


1.222 



TABLE I. Statistical analysis of the critical conductance 
distribution for different boundary conditions (P periodic and 
HW hard wall). Besides the averages of g and lng also the 
standard deviations of these quantities, a g and ai ng are given. 




FIG. 2. As a function of system size the average (ing) is 
plotted for periodic (circles) and hard wall boundary condi- 
tions (squares). 

The average is always over the disorder ensemble. 
Fig.^ shows that the average logarithmic conductance 
still depends slightly on the system size in the regime 
investigated. But the difference between periodic and 
hard wall BCs does not diminish with increasing L, and 
the dependence on L decreases for larger L. Where we 
have used the same system sizes as in Refs. |||9| our 
values for all quantities calculated ((g), (Ing), and the 
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standard deviations of g and ln<?) coincide within one 
percent with the values given in these references. For 
comparison with Q g should be multiplied with a factor 
2, since we consider only one spin direction. Thus, the 
discrepancy between [pi and can indeed be explained 
by the influence of the BCs (see also ficfl). 

Our result has important implications for the scaling 
theory of the metal-insulator transition, since it shows 
that the scaling function (3(g) must depend on the BCs. 
The conductance that enters in this equation has to be 
understood as an average conductance |f|] , and the crit- 
ical conductance is given by (3(g c ) — 0. According to 
our results g c depends on the BCs, g c = 0.413 for peri- 
odic BCs and g c = 0.348 for hard walls at L = 16 (see 
table |), and therefore the [3(g) curves must at least be 
shifted as a function of the BCs. In the next section we 
show by re-examining the weak localization corrections 
to the conductance that also in the metallic regime [3(g) 
depends on the BCs. 



The argument y is defined as 



y 



L 2 



(7) 



Previous analyses in the literature proceeded by approxi- 
mating the sum by an integral []l8f , whereupon all depen- 
dence on the boundary conditions is lost. While this is a 
good approximation for g — > oo, important corrections of 
the order (hig)/g arise for finite g, which we are going to 
derive now, assuming that to this order no further dia- 

§ms beyond the diffuson approximation contribute. In 
it was shown by field theoretical methods combined 
with a renormalization group approach that the diffuson 
approximation gives the leading perturbative contribu- 
tion to the small energy behavior of the spectral correla- 
tion function to order l/g 2 . 

In order to proceed it is convenient to differentiate 
S B c(y)- The derivatives for both PBCs and HW BCs 
can be written with the help of the function 



III. METALLIC REGIME 



(8) 



It is well known that already in metallic regime g> 1 
the quantum interference of diffusing electrons reduces 
the conductance compared to the classical value g = <jL, 
where a is bulk conductivity. The weak localization cor- 
rection 8g is given by a sum over diffusion modes as |lq] 



as 



6g 



-2E 



q 2 L 2 



(3) 



The sum is limited to the diffusive regime where Dq 2 <C 
1/t c . This limitation is taken into account by the expo- 
nential cut-off; r e is the elastic collision time, D = VpT e /3 
denotes the diffusion coefficient and vf the Fermi veloc- 
ity. The sum (|^) depends on the BCs via the quantiza- 
tion condition for the diffusion modes q. For the trans- 
port direction the wave vector is quantized according to 
Qx — n x Tr/L, n x = 1, 2, . . .. Periodic boundary conditions 



in the y-direction imply q y = n y 2n/L, 



±1,±2,. 



and correspondingly for the z-direction. Hard wall BCs 
on the other hand lead to q y = n y Tr/ L, n y = 0,1,2,... 
and q z = n 2 7r/L, n z — 0, 1, 2, . . .. Consequently we have 



Sg 



;S B c(y) 



(4) 



where the index BC stands for a boundary condition and 



Sp(y)= £ 



exp [-7r 2 (n 2 + An y + ^n 2 z )y] 
n 2 , + 4n 2 + An 2 



^>0 

n v ,n z >0 



(5) 
(0) 



d y S P (y) = -n 2 F(y)(2F(4y)) 2 , 
d y S HW (y) = -7T 2 F(y)(l + F(y)) 2 



(9) 
(10) 




The function F(y) is related to the complete elliptic in- 
tegrals K = K{k) and K' = K(k') with k' = s/i - k 2 by 



(11) 



Since we are interested in y <C 1, we need K'/K <C 1 
and therefore k — > 1 (k' -c 1). For small values of k' the 
elliptic integrals behave like 

K = ln^ + 0(k' 2 ),K< = ^ + 0(k' 2 ), (12) 

and we therefore obtain 




(13) 



Inserting this into eqs.(^|), ( jTo| ) and integrating with re- 
spect to y yields 

S P = ^ + \-K\ny-2^ 2 ^+^ y ~a P (14) 
4yfy 8 2 

■\pK 1 1 7T 2 

Shw = -r -= - ^7rhi?/+ -ir 3/2 y/y + — y- a HW . (15) 
A^/y 8 4 8 

whereas, replacing the sum (^) by an integral, one would 
have found 
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s 



a . 



(16) 



where a is an integration constant resulting from the cut- 
off at small q ~ 1/L. Thus, the leading term for small 
V j \/ 7r /y/4, is the same for both boundary conditions. 
The integration constants ctp and a/jw can be evaluated 
numerically, by subtracting from the exactly calculated 
sums the analytical formulae ( |l4|) and (|l5| ) without the 
constants. At the same time this serves as a sensitive 
check for the correctness of these formulae. For small y 
the differences converge to 



a P ~ -6.1509, a w ~ 2.3280 . 



(17) 



We have evaluated the sum numerically down to values 
y = 10~ 6 , where in particular the logarithmic term with 
the prefactors given above could be clearly verified. 
With Eqs.(|i^) and (|l5|) the conductance as a function of 
the dimensionless length L = L/l e takes the form 



g = (cr- A)L - ahiL + b + £>(1/L) 



(18) 



for both periodic and hard wall BCs. The dimensionless 
bulk conductivity a is defined as a — al c h/e 2 , and the 
constant A = \fi/ (2-k 3/2 ) is the same for both BCs. The 
coefficients a and b on the other hand do depend on the 
boundary conditions; their values are given in table 0. 
Note that in the traditional approach the coefficient a 
vanishes! 

Quite surprisingly a < for PBCs, which means that 
the conductance increases even slightly faster than lin- 
early with the system size. This looks as if there was 
anti-localization, but it should be noted that the lead- 
ing behavior due to weak localization is still the usual 
decrease of the (bulk-) conductivity, i.e. the leading term 
is linear in the system size and with the expected neg- 
ative sign. The fact that a < only for PBCs suggests 
a simple physical explanation for the logarithmic term: 
Closing the sample to a double torus by imposing PBCs 
allows for additional paths that interfere constructively 
and lead to enhanced localization for small system sizes 
compared to the HW case. When increasing the system 
size these additional localizing paths quickly stop con- 
tributing and the conductance therefore increases more 
rapidly than what would be expected just from the vol- 
ume part of the weak localization. 

We are now in the position to explore the consequences 
of the BC dependent weak localization corrections for the 
scaling function (3(g). Inserting ([18]) into the definition 



0(g) = 



dhig 
d\nL 



yields 



p{g) = l+-[ahxL-b-a + 0(l/L) 



(19) 



(20) 



It remains to reexpress L by g. To this end we invert 
g(L) from (|l8| ) to order 1/g, 



L 



1 



(g + a In g - a ln(5' - A) + b) , (21) 



a- A 

and insert it in (^0|). We obtain the final result 

(3(g) = 1 - - (b + o(l + ln(er - A)) - alng) + 0(1/ g 2 ) . 
9 

(22) 

It is now obvious that the scaling function does indeed 
depend on the BCs via the coefficients a and b, and 
the dependence arises at order (In g)/g. Furthermore, 
(3(g) depends to order 1/g as well on the material depen- 
dent dimensionless bulk conductivity cf, and is therefore 
non-universal! Again, the non-universality vanishes for 
g — > oo (equivalcntly, on the metallic side of the transi- 
tion: L — > oo), but is important if one is interested in 
(3(g) at finite values of g. Since HW BCs lead to smaller 
values of (3(g) at intermediate values of g than PBCs but 
to a smaller critical conductance, there should be a point 
where the two curves cross, which would imply that in 
that point the change of g with the system size is inde- 
pendent of the BCs. Due to the dependence of (3(g) on 
a, this point is not expected to be universal, though. 

The most interesting question is of course, whether also 
the slope of (3(g) at g = g c is changed by the BCs and/or 
a, as this slope determines the critical exponent v defined 
by i(w) oc Iff - g c \- v according to (3(g) = \(g - g c )/g c . 
This question arises actually already from the depen- 
dence of spectral statistics on the BCs, since the scal- 
ing function can be determined also from purely spectral 
statistics Jl6|,[l7]]. Very recently it has been argued that 
within the same universality class v does at least not de- 
pend on the shape of the sample (2i| . Since the critical 
spectral statistics does depend on the shape of the sam- 
ple much in the same way as on the BCs (indeed, all 
that has been said above about the dependence on the 
BCs translates one to one to a dependence on the shape 
of the sample), one might suspect that v is also inde- 
pendent of the BCs. On the other hand, considering the 
qualitative behavior of the two scaling curves a critical 
exponent independent of the BCs would appear rather as 
coincidence. However, so far it is an open question and 
definitly deserves attention [p6[ . 

With the dependence of the critical conductance dis- 
tribution on the BCs, an experimental test of the correct- 
ness of the one-parameter scaling picture seems within 
reach. Even though an accurate absolute measurement of 
the critical exponent is rather difficult p7|)P8[ , one might 
hope to detect a change with the BCs. To this end it is 
not even necessary to open and close the sample. Rather 
one can investigate the difference between periodic and 
anti-periodic boundary conditions. At least in one direc- 
tion anti-periodic BCs, i.e. a phase factor —1 between 
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two opposite sides of the sample, can be easily produced 
by closing the sample to a ring and introducing half a 
magentic flux quantum (0 = 1/2 in (|l|)). Note that for 
4> = 1/2 the system still belongs to the orthogonal uni- 
versality class, since the Hamiltonian has a real represen- 
tation. This situation has been termed "false time rever- 
sal symmetry breaking" Jl3j |. An experimental search of 
a change of the scaling function in the metallic regime 
upon inclusion of half a flux quantum would be as well a 
most welcome contribution to the long-lasting debate on 
the limits of validity of one-parameter scaling. 

In summary, we have shown that the conductance dis- 
tribution at the Anderson Metal-Insulator transition de- 
pends on the boundary conditions applied in the direc- 
tions transverse to the transport. Furthermore, in the 
metallic regime the dependence of a change of the con- 
ductance with the system size does not depend solely 
on the conductance itself but as well on the boundary 
conditions and the dimensionless bulk conductivity. As 
a consequence the scaling function (3(g) that describes 
the change of conductance when the size of the sample 
is changed is not entirely universal but depends on the 
boundary conditions and the amount of disorder in the 
sample from the metallic regime up to the metal insulator 
transition. 
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a 


b 


PBC 


-5/(2tt) ~ -0.7958 


5(ln3)/(47r) - 2 • 6.1509^ ~ 


-0.8093 


HW 


1/(2tt) ~ 0.1592 


2-2.328/^ - (ln3)/(47r) ~ 


0.3843 



TABLE II. Coefficients in the 1/L expansion of g for periodic boundary conditions (PBC) and hard walls (HW). 
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